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Surface  waves  are  the  roughness  element  of  the  ocean  surface,  the  air-sea  interaction 
processes  are  influenced  by  the  wave  conditions.  The  dynamic  influence  of  surface 
waves  decays  exponentially  with  distance  from  the  air-water  interface.  The  relevant 
length  scale  characterizing  the  decay  rate  is  the  wavelength.  The  parameterization  of 
drag  coefficient  and  surface  roughness  can  be  significantly  improved  by  using  wave¬ 
length  as  the  reference  length  scale  of  atmospheric  measurements.  The  wavelength 
scaling  of  drag  coefficient  and  dynamic  roughness  also  receives  support  from  theo¬ 
retical  studies  of  wind  and  wave  coupling. 
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1.  Introduction 

Drag  coefficient  (CD)  and  surface  roughness  (z0)  are 
important  parameters  for  quantifying  air-sea  momentum 
and  energy  exchanges.  For  a  long  time  the  drag  coeffi¬ 
cient  of  the  ocean  surface  has  been  considered  a  constant 
(e.g.,  Kraus,  1972).  As  the  quality  and  quantity  of  meas¬ 
urements  have  improved,  it  has  become  evident  that  drag 
coefficient  tends  to  increase  with  increasing  wind  speed 
(e.g.,  Garratt,  1977;  Wu,  1980).  Later,  it  was  found  that 
in  addition  to  wind  speed,  stability  stratification,  sea  state 
and  wave  conditions  are  also  important  in  modifying  the 
drag  coefficient  (e.g.,  Donelan,  1982,  1990;  Toba  et  aLy 
1990;  Jones  and  Toba,  2001).  The  study  of  drag  coeffi¬ 
cient  can  also  be  approached  by  investigating  the  surface 
roughness  parameter,  frequently  expressed  in  the 
dimensionless  form  z0*  =  8Zq/u *2,  where  g  is  the  gravita¬ 
tional  acceleration  and  m*  is  the  wind  friction  velocity. 
The  dimensionless  roughness  z0*  is  generally  referred  to 
as  the  Charnock  parameter  and  was  originally  considered 
a  constant  (e.g.,  Charnock,  1955;  Wu,  1980).  Later  on, 
more  extensive  measurements  have  convincingly  illus¬ 
trated  that  z0*  does  not  remain  constant  and  there  seems 
to  be  good  correlation  between  z0+  and  surface  wave  prop¬ 
erties  such  as  wave  height  and  wave  slope  (e.g.,  Donelan, 
1990;  Donelan  et  al.t  1993;  Anctil  and  Donelan,  1996; 
Taylor  and  Yelland,  2001)  or  wave  period  (e.g.,  Toba  et 
a/.,  1990). 

Over  many  decades  of  intensive  research,  our  knowl¬ 
edge  of  the  stability  correction  has  advanced  significantly 
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(e.g.,  Garratt,  1977;  Donelan,  1990;  Geernaert,  1999;  Toba 
et  dl.y  200 1 ).  In  contrast,  a  consistent  parameterization  of 
CD  with  wave  parameters  remains  elusive.  Some  earlier 
works  have  introduced  wave  height  or  wavelength  into 
the  parameterization  of  drag  coefficient  and  surface 
roughness  (e.g.,  Kitaigorodskii,  1973;  Stewart,  1974; 
Donelan,  1990).  Recently,  Makin  and  Kudryavtsev  (1999, 
2002)  and  Makin  (2003)  developed  wind-over-wave  cou¬ 
pling  (WOWC)  theory,  focusing  on  the  dynamics  of 
microscale  air-sea  interactions.  The  surface  wind  stress 
has  been  shown  to  be  significantly  influenced  by  the  sur¬ 
face  wave  dynamics  through  the  form  drag,  which  is  the 
correlation  of  wave-induced  pressure  field  and  the  sur¬ 
face  wave  slope.  The  form  drag  is  further  divided  into  a 
non-separated  sheltering  component  and  the  stress  due  to 
flow  separation  from  breaking  of  both  small-scale  and 
dominant-scale  waves.  Numerical  computations  have 
been  carried  out  to  demonstrate  the  relative  importance 
of  the  different  surface  stress  components  and  their  de¬ 
pendence  on  wind  speed,  wave  age,  finite  water  depth 
and  other  atmospheric  and  wave  parameters.  In  the 
WOWC  analysis,  the  wavelength  at  the  spectral  peak  is  a 
clear  scaling  length  of  air-sea  interaction  processes. 

Almost  all  the  earlier  analyses  of  CD  and  z0  use  the 
reference  wind  speed  at  10  m  elevation  (£/10).  Although 
the  employment  of  Ul0  in  the  analyses  provides  a  con¬ 
sistent  reference  level  of  wind  measurements  (as  com¬ 
pared  to  earlier  reports  using  “mast  height”  or  “anemom¬ 
eter  height”),  the  dynamical  significance  of  the  10-m  el¬ 
evation  in  the  marine  boundary  layer  is  not  clear.  From  a 
heuristic  point  of  view,  surface  waves  are  the  ocean  sur¬ 
face  roughness,  the  air-sea  interaction  processes  are  in¬ 
fluenced  by  the  wave  conditions.  Because  the  influence 
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of  surface  waves  decays  exponentially  with  the  wave¬ 
length  serving  as  the  vertical  length  scale,  the  dynami¬ 
cally  meaningful  reference  elevation  should  be  the  char¬ 
acteristic  wavelength  of  the  peak  component  of  surface 
wave  spectrum,  A.  As  described  in  the  last  paragraph,  this 
common  sense  vertical  length  scale  also  receives  theo¬ 
retical  support  from  the  study  of  wind  over  waves  cou¬ 
pling.  This  paper  presents  an  analysis  with  Um  serving 
as  the  reference  wind  speed,  where  U^2  is  the  wind  speed 
at  the  elevation  equal  to  one-half  of  the  surface  wave¬ 
length.  The  general  conclusions  of  this  paper  are  not 
changed  by  using  different  fractions  of  the  wavelength 
as  the  vertical  length  scale  of  wind  speed  reference;  more 
discussions  on  this  point  will  be  presented  later.  The  func¬ 
tional  expressions  of  the  drag  coefficient  referenced  to 
U XJ2  (expressed  as  C^2  to  distinguish  it  from  C10,  the  drag 
coefficient  referenced  to  Ul0)  and  surface  roughness  and 
their  dependence  on  the  wind  and  wave  parameters  are 
described.  The  dimensionless  roughness  in  the  form  of 
kpzQ  falls  out  naturally  when  wavelength  is  used  as  the 
reference  length  scale,  where  kp  =  2^/A  is  the  wavenumber 
at  the  spectral  peak.  The  relation  between  kpzQ  and  z0* 
can  be  easily  derived. 

The  detail  of  the  analysis  with  Um  as  the  reference 
wind  speed  and  the  parameterization  of  surface  rough¬ 
ness  in  terms  of  kpz0  and  z0 ♦  is  given  in  Section  2.  Field 
data  with  wind  stress  measurements  under  conditions 
dominated  by  wind  seas  are  assembled  to  verify  the  re¬ 
sults  obtained  from  the  analysis  (Section  3).  A  compari¬ 
son  of  and  C10  illustrates  the  superiority  of  the  former 
for  grouping  together  datasets  from  different  sources. 
Similarly,  the  dimensionless  roughness  represented  by 
kpz0  consolidates  data  better  than  z0*.  Section  4  presents 
discussions  on  the  issue  of  possible  spurious  correlation 
when  scaling  with  w*  (Smith  et  at. ,  1992)  in  the 
parameterization  function  of  the  drag  coefficient  and  sur¬ 
face  roughness.  The  analysis  indicates  that  spurious  cor¬ 
relation  does  not  pose  a  problem  in  the  datasets  used  in 
the  present  study.  The  summary  and  conclusions  are  given 
in  Section  5. 

2.  Analysis 

In  this  paper,  the  discussion  is  limited  to  neutral 
stratification  or  to  the  situation  that  the  stability  correc¬ 
tion  has  been  performed.  The  vertical  wind  speed  distri¬ 
bution  is  described  by  the  logarithmic  wind  profile, 

y(z)  =  — In— ,  (1) 

K  Zq 

where  U  is  the  equivalent  wind  speed  at  neutral  condi¬ 
tion,  k*  is  the  von  Karman  constant  (0.4),  and  z  is  the  ver¬ 
tical  elevation  measured  from  the  mean  water  surface. 


Taking  the  reference  length  scale  for  wind  speed  meas¬ 
urements  at  the  level  equal  to  one-half  of  the  wavelength, 
the  drag  coefficient  C^2  =  u^lU^ 2  can  be  written  as 


"A/2 


i 


(2) 


The  approach  of  using  kp  for  the  characteristic  wave  com¬ 
ponent  is  expected  to  work  well  in  wind-sea  dominated 
wave  conditions.  The  influence  of  swell  on  the  drag  co¬ 
efficient  and  surface  roughness  requires  further  investi¬ 
gation. 

Equation  (2)  suggests  that  a  natural  expression  of 
the  dimensionless  surface  roughness  is  kpz0 ,  and  from  (2) 

kpz0=Jtexp(-KCfii).  (3) 


Using  the  dispersion  relation  of  surface  waves. 


£w;;=sVanh*A  (4) 

Mo  can  appear  in  several  different  forms 


Mo  =  z°*[ ^ j  tanh  kPh  =  z°*  ~£~ ] 


C^/2  tanh  kph 


Zp.cot  z0,co?Cx/2 
tanh  kph  tanh  kph 


(5) 


where  £0*.  =  uta)p/g  and  ft),  =  UU2(oplg,  both  can  be  inter¬ 
preted  as  the  dimensionless  frequency  of  the  air-sea  cou¬ 
pled  system;  cop  and  Cp  are  the  angular  frequency  and 
phase  speed,  respectively,  of  the  peak  spectral  compo¬ 
nent.  For  deep  water,  the  dimensionless  frequency  is  iden¬ 
tical  to  the  inverse  wave  age  uJCp  or  Um!Cp. 

The  analysis  of  the  logarithmic  wind  profile  leads  to 
one  equation  (Eq.  (2)  or  (3))  connecting  two 
dimensionless  unknowns  ( Cm  and  kpz0).  To  make  further 
progress,  either  another  wind-wave  coupling  function  or 
assumption  needs  to  be  introduced  or  an  empirical  rela¬ 
tion  of  Cm  or  kpz o  needs  to  be  established.  The  next  sec¬ 
tion  presents  an  analysis  of  results  from  several  field  ex¬ 
periments  (Donelan,  1979;  Merzi  and  Graf,  1985;  Anctil 
andDonelan,  1996;  Janssen,  1997).  The  selected  datasets 
are  characterized  by  wind-sea  dominated  wave  conditions 
and  the  tabulated  data  list  sufficient  information  on  wave 
properties  and  water  depth  to  derive  the  wavelength.  The 
investigation  shows  that  among  the  dependence  of  CV2, 
kpZp  and  z0,  on  uJCp,  UXn/Cp,  a )„  or  g>„  the 
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Fig.  1.  (a)  Drag  coefficient  Cm((o„)  measured  in  wind-sea  dominated  wave  conditions.  The  numbers  within  the  square  brackets 
in  the  legends  are  the  maximum  and  minimum  wavelengths  in  each  individual  dataset;  the  curve  represents  the  parameterization 
equation  (6).  (b)  Dimensionless  roughness  in  terms  of  kpz0(co„),  the  curve  represents  the  parameterization  equation  (7). 
(c)  Dimensionless  roughness  in  terms  of  z0*(fiO,  the  curve  represents  the  parameterization  equation  (8).  (c)  also  shows  the 
curve  from  least  square  power-law  fitting  of  z0.(Cyw*).  The  functions  given  by  Maat  et  al.  (1991)  and  Smith  et  al.  (1992)  are 
also  illustrated  (Section  4)  for  comparison. 


parameterization  of  yields  a  simple  power-law 

function, 

^A/2  =  Accq*z  ,  (6) 

where  Ac  =  1.22  x  10~2  and  ac  =  0.704.  More  details  of 
the  analysis  of  the  field  data  are  presented  in  the  next 
section.  The  dimensionless  roughness  becomes 

kpz0=7t  exp|^r(4ot#  )~°'5  J  (7) 

The  dimensionless  roughness  expressed  as  the  Charnock 
coefficient  can  be  obtained  from  substitution  of  (5)  in  (7), 

z0 *  =  n  tanh  kh(0~}  exp^-K^GJ?'  )~°’5  j.  (8) 

Other  expressions  of  the  dimensionless  roughness  in  terms 
of  uJCp  or  Um!Cp  can  also  be  constructed  and  will  not 
be  further  elaborated  here. 


3,  Comparison  with  Field  Measurements 

In  order  to  verify  the  analysis  presented  in  the  last 
section,  tabulated  field  data  reported  by  Donelan  (1979), 
Merzi  and  Graf  (1985),  Anctil  and  Donelan  (1996)  and 
Janssen  (1997)  have  been  assembled.  All  these  datasets 
share  several  common  attributes,  including  data  quality 
control  screening  for  steady  wind  events,  fetch-limited 
wave  conditions,  direct  wind  stress  measurements  using 
either  the  profile  method  or  the  eddy  correlation  tech¬ 
nique,  and  they  contain  sufficient  information  to  derive 
the  peak  wavelength.  The  listed  wind  velocity  is  corrected 
for  the  stability  influences.  Moreover,  most  of  the  data 
have  been  collected  in  finite-depth  waters;  therefore,  a 
proper  accounting  of  the  depth  effect  as  described  in  (5) 
or  (8)  should  be  applied  when  the  Charnock  parameter  is 
used  for  the  dimensionless  surface  roughness.  However, 
the  data  scatter  is  rather  large,  the  depth  effect  is  masked 
and  the  scatter  plots  of  z0*(co **)  and  z0*/tanh(*/0(G>*,)  are 
not  very  different  (not  shown).  A  brief  description  of  these 
experiments  is  given  in  Appendix  A. 

Figure  1(a)  displays  the  dependence  of  Cy2  on 
processed  from  the  assembled  data.  Although  the  data 
scatter  is  large,  the  result  is  a  substantial  improvement 
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over  the  expressions  of  C10.  Figure  A1  in  Appendix  A 
shows  several  presentations  of  C10,  including  (a)  CI0(t/10), 
(b)  C10(ft)* *),  and  (c)  Cl0(U]0ICp).  The  C10  representations 
result  in  sets  of  curves  and  clearly  require  a  secondary 
parameter.  Comparing  the  results  presented  in  Figs.  1  and 
A 1,  it  is  obvious  that  Cy2(cot+)  does  indeed  provide  a  bet¬ 
ter  parameterization  than  C,0  regarding  the  surface  wave 
effects  on  the  drag  coefficient.  Polynomial  fitting  of  a 
power  law  function  to  the  combined  data  presented  in  Fig. 
1(a)  yields  the  empirical  function  of  Cy2(CQt*)  given  in 
(6).  (Different  parameterizations  of  Cm  with  uJCp  and 
other  combinations  of  wind  and  wave  parameters  have 
also  been  investigated.  For  example,  the  regression  con¬ 
stant  and  exponent  for  uJCp  parameterization  are  1.09  x 
10’2  and  0.696,  respectively.  The  root  mean  square  (rms) 
difference  between  the  regression  curve  and  the  meas¬ 
urements  is  2.78  x  lO"4,  which  is  about  20  percent  larger 
than  Cy2((t)„)  parameterization  with  an  rms  difference 
of  2.26  x  10"4.)  The  issue  of  the  fraction  of  wavelength 
as  the  vertical  reference  level  is  also  studied.  Examples 
of  using  Uy 4,  Uy6>  Uy 8  as  reference  wind  speeds  are  given 
in  Appendix  B,  showing  that  the  wavelength  scaling  of 
the  drag  coefficient  is  not  influenced  by  the  choice  of  the 
proportionality  constant. 

The  dimensionless  roughness  kpz0(co„)  is  shown  in 
Fig.  1  (b).  The  measured  data  again  display  a  smooth  vari¬ 
ation  with  ft)**  and  measurements  from  different  sources 
collapse  very  well  following  the  proposed 
parameterization  function  (7),  given  as  the  solid  curve  in 
the  figure.  In  terms  of  the  Charnock  parameter,  z0*  (Fig. 
1(c)),  the  dimensionless  roughness  increases  with  ft)**  in 
the  parameter  range  of  the  available  data.  The 
parameterization  function  (8)  predicts  that  the  depend¬ 
ence  of  z0*(ft)**)  reverses  to  a  decreasing  trend  at  ft)**  > 
-0.25  and  z0*  approaches  ft)**’2  asymptotically  as  ft)**  — > 
00 •  Field  data  at  very  high  ft)**  are  difficult  to  obtain  so 
such  non-monotonic  behavior  cannot  be  positively  con¬ 
firmed  from  the  field  measurements  assembled  here.  The 
inverse  relationship  of  z0*(ft)**)  at  high  ft)**  range  has  been 
observed  in  laboratory  measurements  (e.g.,  Toba  et  al. , 
1990). 

4.  Discussions 

From  the  point  of  view  of  dimensional  analysis,  (6) 
can  be  interpreted  as  choosing  peak  wavelength  and  phase 
velocity  of  surface  waves  as  the  length  and  velocity  scales 
for  organizing  the  observations.  As  discussed  earlier,  these 
scaling  factors  are  logical  choices  because  the  roughness 
of  the  ocean  surface  is  primarily  composed  of  surface 
waves.  Moreover,  the  study  of  wind  over  wave  coupling 
supports  the  suggestion  of  peak  wavelength  as  the  verti¬ 
cal  scaling  length  for  parameterizations  (Makin  and 
Kudryavtsev,  1999,  2002;  Makin,  2003). 

Smith  et  al.  (1992)  caution  about  the  possible  haz¬ 


ards  of  self-correlation  in  scaling  with  «*  because  the 
exponential  dependence  of  z0  on  m*  (and  that  wind  speeds 
at  different  elevations — e.g.,  from  instrument  height  to 
Z  =  10  m  or  A/2 — are  subsequently  computed  from  z0  and 
w*  using  the  function  describing  the  wind  speed  profile). 
They  provide  two  test  criteria  for  determining  whether 
the  spurious  correlation  is  negligible  in  a  dataset  (p.  134, 
ibid.)  with  z0*  -  (C/u*)a, 

var(lnZ)  «  var(InF)*  and  var(lnW)  «  var(lnA'),  (9) 

where  Z=ua,Y=  (gz/u2)txp(-KU/ut)u*a,  W  =  w*,  X  =  Cp 
and  the  function  var  is  the  variance  of  the  argument  in 
parentheses.  They  further  comment  that  the  second  con¬ 
dition  is  usually  not  met  in  most  field  data  (in  the  exam¬ 
ple  cited  by  them,  the  four  variances  are  var(InZ)  =  0. 146, 
var(lnT)  =  0.486,  var(lnW0  =  0.146,  and  var(lnX)  = 
0.0057),  thus  “the  variability  of  lnw*  is  not  smaller  than 
that  of  InCj,,  variation  of  \n(Cp/u *)  is  due  more  to  w*  than 
to  Cpy  and  self-correlation  has  influenced  the  fit  to  obtain 
(the  parameterization  equation).”  They  further  comment 
that  “This  does  not  disprove  (the  parameterized  equation) 
...”  and  that  “As  far  as  we  know,  other  similar  data  sets 
would  lead  to  the  same  impasse  if  tested  for  spurious¬ 
ness.” 

The  four  datasets  assembled  in  this  paper  cover  a 
considerably  wider  range  of  wind  sea  conditions.  The  least 
square  fitting  of  the  combined  data  yields  z0*  =  1 .50(0^/ 
«*)-1*36.  The  fitted  power-law  function  is  shown  in  Fig. 
1(c)  in  terms  of  z0*(ft)**),  where  Cp/u *  =  ft)**’1  for  deep 
water  condition.  As  noted  earlier,  the  large  data  scatter 
has  masked  the  depth  effect  in  these  datasets.  For  com¬ 
parison,  the  functions  presented  by  Maat  et  al.  (1991)  and 
Smith  et  al.  (1992),  z0*  =  0.8 (C/w*)’1  and  z0*  =  0.48(0/ 
h*)’1,  respectively,  are  also  illustrated.  The  computed 
variances  are  var(lnZ)  =  0.278,  var(lnK)  =  0.687, 
var(InVf)  =  0.150,  and  var(lnX)  =  0.202.  The  problem  of 
spurious  self-correlation  is  not  expected  to  be  serious  in 
the  present  analysis. 

5.  Summary  and  Conclusions 

The  dynamic  influence  of  surface  waves  decays 
exponentially  with  distance  from  the  air-water  interface 
and  the  decay  rate  is  inversely  proportional  to  the  char¬ 
acteristic  wavelength.  These  dynamic  properties  of  sur¬ 
face  waves  suggest  that  the  reference  wind  speed  for  at¬ 
mospheric  measurements  on  air-sea  coupling  parameters 
should  be  the  wind  speed  at  an  elevation  proportional  to 
the  wavelength.  The  wavelength  scaling  also  receives 
strong  support  from  wind  over  wave  coupling  theory  (e.g., 
Makin  and  Kudryavtsev,  1999,  2002;  Makin,  2003).  The 
analytical  expressions  of  the  drag  coefficient  and  surface 
roughness  with  Uy2  as  reference  wind  speed  are  estab¬ 
lished  from  the  logarithmic  wind  profile  (Section  2)  for 
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Fig.  AI.  Drag  coefficient  in  terms  of  C10  presented  as  (a)  C10((/l0),  (b)  Cl0(a>..),  and  (c)  CX0(Ui{JCp).  The  data  scatter  in  these 
expressions  is  considerably  worse  than  the  <*>..)  expression  shown  in  Fig.  1(a)  and  suggests  the  necessity  of  a  secondary 
parameter  for  C!0. 


stratification-corrected  conditions  and  wind-sea  domi¬ 
nated  wave  fields.  The  analysis  suggests  that  a  natural 
expression  of  the  dimensionless  surface  roughness  is  kpz0 . 
Other  expressions  of  the  nondimensional  surface  rough¬ 
ness,  such  as  the  Charnock  parameter,  can  be  obtained  by 
the  substitution  of  kp  with  cop  using  the  dispersion  rela¬ 
tion  of  surface  gravity  waves. 

Several  field  experiments  with  wind-sea  dominated 
wave  fields  and  tabulated  data  that  provide  sufficient  in¬ 
formation  to  derive  wavelengths  have  been  analyzed.  The 
results  show  that  the  measurements  of  the  drag  coeffi¬ 
cient  from  different  experiments  can  be  consolidated  in  a 
more  organized  manner  when  it  is  presented  in  terms  of 
Cm  than  in  terms  of  C10  (compare  Figs.  1(a)  and  Al). 
With  the  consideration  that  wavelength  and  wave  phase 
speed  are  the  proper  length  and  velocity  scales  of  the 
ocean  surface  roughness,  the  parameterization  of  Cm(co**) 
is  given  by  an  empirical  power-law  function  (6).  The 
dimensionless  roughness  can  then  be  represented  by  (7) 
for  kpz0(G )**)  or  by  (8)  for  z0*(ft)**).  Surface  roughness  in 
terms  of  kpz0  increases  monotonically  with  ft)**.  When 
expressed  in  terms  of  z0*  the  roughness  increases  with 
ft)**  for  ft)**  <  -0.25  and  decreases  with  ft)**  for  ft)**  > 
-0.25.  The  range  of  ft)**  in  the  field  data  is  usually  much 
less  than  0.25  but  laboratory  data  extend  to  a  much  higher 
ft)**  range  (e.g.,  Toba  et  al 1990). 
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Appendix  A.  Descriptions  of  the  Assembled  Datasets 

Donelan  (1979)  acquires  wind  and  wave  measure¬ 
ments  from  a  fixed  tower  at  the  western  side  of  Lake 
Ontario.  The  tower  is  located  1 100  m  from  the  coast,  and 
the  local  water  depth  is  12  m.  An  extensive  suite  of  in¬ 
struments  are  mounted  on  the  tower.  Data  relevant  to  this 
study  include  wind  velocity,  direct  wind  stress  measure¬ 
ment  using  the  eddy  correlation  method,  wave  variance 
and  peak  frequency  of  the  wave  spectrum.  The  wind  con¬ 
ditions  selected  for  his  analysis  are  steady  offshore  blow¬ 
ing  events  with  directions  of  both  winds  and  waves  stay¬ 
ing  within  25  degrees  to  the  beach  normal,  and  that  the 
peak  frequency  of  a  companion  waverider  buoy  deployed 
further  offshore  to  be  less  than  3.14  rad/s.  All  together, 
six  cases  are  reported.  The  maximum  and  minimum  wave¬ 
lengths  in  this  dataset  are  5.7  and  3.3  m. 

Merzi  and  Graf  (1985)  conduct  wind  stress  meas¬ 
urements  in  Lake  Geneva.  The  measurement  station  is  on 
the  northern  part  of  the  lake.  The  local  water  depth  is  3 
m.  The  records  selected  for  analysis  are  based  on  steady 
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Fig.  BI.  Comparison  of  (a)  C (b)  C^4(comm),  (c)  C^6((Omm)y  and  (d)  Cm{(0 **).  The  wavelength  scaling  for  the  drag  coeffi¬ 
cient  is  not  qualitatively  influenced  by  the  choice  of  the  proportionality  constant. 


southwesterly  events  and  well-behaved  wind  profiles  (60 
cases).  The  friction  velocity  is  derived  from  the  profile 
method.  Thermal  stratification  is  corrected.  The  friction 
velocity  used  in  this  analysis  is  calculated  from  their  tabu¬ 
lated  neutral  drag  coefficient  and  Ul0.  The  maximum  and 
minimum  wavelengths  in  this  dataset  are  19.0  and  8.6  m. 

Anctil  and  Donelan  (1996)  report  momentum  flux 
observations  from  four  towers  at  different  depths  along  a 
shore-normal  line  at  the  west  end  of  Lake  Ontario, 
Canada.  The  nominal  depths  of  the  towers  are  2,  4,  8,  and 
12  m,  the  wind  stress  is  measured  by  the  eddy  correlation 
method  at  anemometer  heights  of  6.2,  7.8,  7.4  and  12  m 
for  the  four  towers.  The  drag  coefficients  and  wind  speeds 
at  both  anemometer  heights  and  10  m  reference  elevation 
are  tabulated  along  with  the  stability  parameter,  instanta¬ 
neous  water  depth,  root  mean  square  wave  height,  peak 
wave  period,  peak  phase  speed,  wave  angle  and  steep¬ 
ness.  All  together,  18  cases  are  reported.  The  maximum 
and  minimum  wavelengths  in  this  dataset  are  75.8  and 
28.6  m. 

Janssen  (1997)  carries  out  a  statistical  error  analysis 
of  the  HEXMAX  data  to  investigate  the  sea-state  influ¬ 
ence  on  wind-stress.  The  HEXMAX  data  are  collected 
from  a  research  platform  at  a  nominal  depth  of  18  m.  The 
wind  stress  is  calculated  by  the  eddy  correlation  method. 
The  dataset  (58  cases)  has  been  carefully  screened  to  re¬ 
move  unsteady  events  and  swell-contaminated  cases.  The 
wind  stress,  wind  speed  at  10  m  elevation,  peak  phase 
speed  and  significant  wave  height  are  tabulated.  The  char¬ 
acteristic  wavelength  needed  for  the  present  analysis  is 


derived  by  the  dispersion  relation  using  the  tabulated 
phase  speed  and  the  nominal  water  depth.  The  maximum 
and  minimum  wavelengths  in  this  dataset  are  101.6  and 
53.4  m. 

Figure  A 1  shows  several  presentations  of  C10,  includ¬ 
ing  (a)  Cl0(t/I0),  (b)  C10(tu*,),  and  (c)  C[0(UxotCp).  The 
dimensionally  consistent  expressions  (Figs.  A 1  (b)  and  (c)) 
reveal  the  stratification  of  Cl0  with  wavelength.  The  wave¬ 
length  ranges  of  the  four  datasets  are  (101.6  m,  53.4  m) 
in  Janssen  (1997),  (75.8  m,  28.6  m)  in  Anctil  and  Donelan 
(1996),  (19.0  m,  8.6  m)  in  Merzi  and  Graf  (1985),  and 
(5.7  m,  3.3  m)  in  Donelan  (1979).  For  dimensionally  in¬ 
consistent  expression  such  as  C10(£/10)  shown  in  Fig. 
A  1(a),  such  a  wave  effect  is  still  discernable.  The  results 
shown  in  Fig.  A1  suggest  that  the  drag  coefficient  in  terms 
of  C10  cannot  be  expressed  as  a  single-parameter  func¬ 
tion.  Two-parameter  functions,  such  as  C10(f/l0,  uJCp) 
and  C10(£/10,  Ul0/Cp)  have  been  reported. 

Appendix  B.  A  Comparison  of  C^2  and  Cxjn 

The  drag  coefficient,  being  the  square  of  the  ratio 
between  the  wind  friction  velocity  and  a  reference  wind 
speed,  is  quantitatively  affected  by  the  choice  of  the  ref¬ 
erence  wind  speed.  It  is  demonstrated  in  Figs.  1  and  A1 
that  the  same  group  of  data  can  be  better  organized  when 
presented  in  terms  of  C^2  in  comparison  with  C,0.  If  an 
elevation  other  than  one-half  of  the  wavelength  is  used 
for  the  vertical  length  scale,  the  equation  of  the  drag  co¬ 
efficient  can  be  written  as 
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c”“=1'n^J  •  <BI) 

where  kp  =  Inti.  The  drag  coefficient  Cm{co..),  Cw(0)„), 
an<^  Cys((Ott)  of  the  assembled  data  are  dis¬ 
played  in  Fig.  Bl.  As  expected,  the  collapsing  of  meas¬ 
urements  is  not  influenced  by  the  choice  of  the  propor¬ 
tionality  constant  for  wavelength  scaling. 
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